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Background: The deviation between different model calculations that may occur when one goes toward regions 
where the masses are unknown is getting increased attention. This is related to the uncertainties of the different 
models which may have not been fully understood. 

Purpose: To explore in detail the effect of the isospin dependence of the spin-orbital force in the Woods-Saxon 
potential on global binding energy and deformation calculations. 

Method: The microscopic energies and nuclear deformations of about 1850 even-even nuclei are calculated sys¬ 
tematically within the macroscopic-microscopic framework using three Woods-Saxon parameterizations, with dif¬ 
ferent isospin dependences, which were constructed mainly for nuclear spectroscopy calculations. Calculations are 
performed in the deformation space (/92,7, Pa)- Both the monopole and doubly stretched quadrupole interactions 
are considered for the pairing channel. 

Results: The ground state deformations obtained by the three calculations are quite similar to each other. Large 
differences are seen mainly in neutron-rich nuclei and in superheavy nuclei. Systematic calculations on the shape- 
coexisting second minima are also presented. As for the microscopic energies of the ground states, the results are 
also very close to each other. Only in a few cases the difference is larger than 2 MeV. The total binding energy is 
estimated by adding the macroscopic energy provided by the usual liquid drop model with its parameters fitted 
through the least square root and minimax criteria. Calculations are also compared with the results of other 
macroscopic-microscopic mass models. 

Conclusions: All the three calculations give similar values for the deformations, microscopic energies and binding 
energies of most nuclei. One may expect to have a better understanding of the isospin dependence of the spin- 
orbital force with more data on proton- and neutron-rich nuclei. 

PACS numbers: 21.10.Dr, 21.30.Fe, 21.60.Jz, 24.10.Cn 


I. INTRODUCTION 

Nuclear physics is an emergent phenomenon where reg¬ 
ular and simple patterns can be created by the com¬ 
plicated interplay among its constitutes: protons and 
neutrons [T|. The understanding of its emergent behav¬ 
ior progresses by systematic experimental observations 
and the construction of models to interpret them. New 
physics has been continuously revealed through the un¬ 
derstanding of the global smooth behaviors of nuclear 
structure as well as its local fluctuations. In particular, 
studies on the nuclear mass and other ground state prop¬ 
erties reveal strikingly systematic behaviors including the 
nuclear liquid-like property, shell structure, pairing cor¬ 
relation and superfluidity as well as the nuclear deforma¬ 
tion. It is not surprising that the so called macroscopic- 
microscopic (mac-mic) model , which is in line with 
above picture, has been extremely successful in describ¬ 
ing the nuclear ground state properties and spectroscopy. 
A reliable and precise theoretical prediction is also of par¬ 
ticular importance for the estimation of the masses of 
experimentally unknown nuclei far from stability, partic¬ 
ular those along the astrophysical r-process path. The 
mac-mic model and its derivations are still extensively 
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applied nowadays [§-[^. With phenomenological cor¬ 
rections included, microscopic nuclear density functional 
calculations can also give a comparable description of the 
binding energy pHI ]. Beyond mean-field corrections are 
also considered in recent global calculations in,HI. Be¬ 
sides these global calculations, local mass formulas and 
the nuclear shell model can give even more precise de¬ 
scription for a limited number of nuclei (see, e.g., Refs. 
[281 - 1^ and references therein). 

One perspective that is getting increased attention is 
the deviation between different effective theories that 
may occur when one goes toward regions where the 
masses are unknown. This is related to the uncertain¬ 
ties of the different models which may have not been 
thoroughly understood [3l|, For example, a num¬ 

ber of parameters of the Skyrme force still can not be 
fully determined by fitting to available experimental data 
and show large uncertainties One encounters 

the same problem when determining the isospin depen¬ 
dence of the spin-orbit (SO) force of the phenomenolog¬ 
ical Woods-Saxon (WS) potential, which can have sig¬ 
nificant effects on the evolution of the shell structure in 
light neutron-rich nuclei [s^ . Such uncertainties can be 
related to the description of the single-particle structure 
in both self-consistent mean field and phenomenological 
approaches [H, The influence of the SO coupling of 
the Skyrme force on global binding energy calculations 
was recently studied in Ref. [^ . 
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The motivation of this work is to calculate system¬ 
atically the microscopic energies and nuclear deforma¬ 
tions within the mac-mic framework with different WS 
parameterizations which were constructed primarily for 
nuclear spectroscopy calculations. In particular, we are 
interested to explore in detail the performance of the WS 
parameterization of Ref. [s^ in heavy nuclei and to see 
whether the isospin dependence of the SO force has any 
global influence on binding energy and deformation. Cal¬ 
culations will also be compared with the results of other 
mac-mic mass models. 

The paper is organized as follows: In Sec. II, we briefly 
introduce the mac-mic approach and the empirical WS 
single-particle potential. It is followed by the detailed 
comparisons on the calculated nuclear deformations and 
microscopic energies in Sec. III. The total binding ener¬ 
gies are studied in Sec IV where the parameters of the 
liquid drop model are fitted to experimental binding en¬ 
ergies. A short summary is given in Sec. V. 


II. THE THEORETICAL FRAMEWORK: 
MICROSCOPIC ENERGY AND THE NUCLEAR 
DEFORMATION 


which correspond to the central potential, the SO poten¬ 
tial and the Coulomb potential, respectively. The surface 
is defined as 


R = C{p)Ro 


A,^ 




( 4 ) 


where C{I3) is the volume conservation factor. The Lund 
convention has been used to transform the coefficients 
to nuclear deformations in terms of /3 and 7 . In the 
present work, a non-axially symmetric shape with coeffi¬ 
cients up to / 34 k is considered but only even multiples are 
included. That is, the potential energy surface calcula¬ 
tion is performed in the deformation space (/ 32 , 7 ,/ 34 ) 

The ground state deformation values are taken as those 
that correspond to the minimum in the total energy. The 
influence of the octupole and higher-order deformations 
is considered later in a systematic calculation within a 
axially symmetric deformation space. 

A variety of parameterizations of the WS potential ex¬ 
ists (see, e.g.. Refs. (dSl - E^ . Table II in Ref. and 
Table I in Ref. [1^). In the “standard” one the 

strengths of the central and SO potentials are given as 


Within the mac-mic framework 0 , the total energy of 
a nucleus can be written as the sum of a macroscopic and 
a microscopic terms. The macroscopic term describes the 
bulk properties of the nucleus. It is usually approximated 
by the standard liquid drop model or its revised versions. 
The microscopic term, which may show large fluctuations 
with changing deformation and particle number, has its 
origins in the quantum shell effects. It describes the sin¬ 
gle particle properties of the nucleons near the Fermi sur¬ 
face. It usually consists of the shell and pairing correction 
terms, which are evaluated in an average potential well. 
Thus, for a given nucleus with Z protons and N neutrons 
at the full set of deformation parameters /3, the total en¬ 
ergy can be written as 

E{N,Z,P) = E^ac{N,Z) + Edef{N,Z,P) 

+ E,hell{N,Z,l3)+Eparr{N,Z,^) ( 1 ) 

The deformation correction energy is written as im 
Eaef = [Bsm - + [B,{P) - 1]A(o) (2) 


where Bg and B^ are functions of the shape of the nucleus 
only and are equal to 1 when the nucleus is spherical. 
E^^'^ and E^^'^ refer to the spherical surface energy and 
the spherical Coulomb energy, respectively. 

In this paper, the single-particle level is derived from 
a non-axial deformed WS potential [ 1 , 0 , of the 

form: 


V{r,l3) 


_ V _ 

1 -I- exp[{r — R)/a] 

^ 1 -h exp[{r - R)/aso] ^ 

+ 2 ^Coul 


( 3 ) 


V = -Vo{l + ^t-Td), 


( 5 ) 


and 


Vso = -XVo{l + ^^t-Td), ( 6 ) 

respectively, where we have replaced the original N — Z 
term with 4t • to get a consistent description of both 
protons and neutron orbitals, t and denote the isospin 
quantum numbers of the last nucleon and of the daughter 
nucleus, respectively. The total isospin of the system is 
T = I-I-Tai-i. It is 4t-T^_i = —3 for the T — 0 ground 
state of a iV = Z nucleus and 

4t • Ta -1 = N — Z — 1 for neutron orbits 

= —{N — Z + 3) for proton orbits (7) 

in N > Z nuclei with T = {N — Z)/2 [d^. In Ref. (d^ . 
the isospin-dependent terms are parameterized as 

33 

n = nso = -4, ( 8 ) 

51 

where the SO potential depth is assumed to have the 
same isospin dependence as that of the central potential. 
This assumption is rather commonly used (d^ . The typ¬ 
ical strength of k is in the range —0.6 ^ —0.9. kso is 
assumed to be zero in Ref. [d^. In Ref. [1^, an uncon¬ 
ventional assumption is taken as 


KSO = -K (9) 

in order to explain the shell evolution in light neutron- 
rich nuclei. We have shown that above choice of kso is 
essential for the description of the disappearence of the 
N = 8 , 20 shell closures and the emergence of the N = 6, 
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14, 16, 32 as well as 34 new subshells in neutron-rich nu¬ 
clei. A number of selected spherical as well as deformed 
nuclei are evaluated in Ref. but its global perfor¬ 
mance has not yet been explored. In particular, our cal¬ 
culations reproduce well the intruder configurations and 
ground state deformations of neutron-rich nuclei around 
= 20. A similar island of inversion is expected in 
neutron-rich nuclei around N = 40. 

The shell correction energy has been calculated by the 
traditional Strutinsky methods [i^. A correction poly¬ 
nomial of the order p = 6 and the smoothing parameter 
7 = 1.2hu}Q{fkjJo = 41A“^/^) has been used in this pa¬ 
per. The pairing energy is calculated considering both 
the monopole and doubly stretched quadrupole interac¬ 
tion: 






where 


( 10 ) 


(V) ^ f _ 

1 {alQ.m 


X — 0, p = 0, 

A = 2,/r = 0,1,2. 


( 11 ) 


Here, a{a) denotes the states of signature r = —i {r = 
—i). Using the double-stretched operator (Q" = r"‘^Y 2 ^), 
the generators of quadrupole pairing interaction in Eqs. 
( 11 ) have the form: 


Qo — Q20J 

Q^ = ^(Q'V + Q"_^), = (12) 


The monopole paring strength, Gqo, is determined by 
the average gap method [b^, while the quadrupole pair¬ 
ing strengths, G 2 i_i, are obtained by restoring the Galilean 
invariance of the system under A-pole collective shape 
oscillations [5l|. In order to avoid the spurious pairing 
phase transition and particle number fluctuation encoun¬ 
tered in the BCS calculatiom the pairing is treated by the 
Lipkin-Nogami approach ^ in which the particle num¬ 
ber is conser ved approximately. For further details, see 
Refs. 



FIG. 1. (Color online) Systematic calculations on the 
quadrupole deformation ^2 of even-even nuclei by using the 
Woods-Saxon potential with WSl from Ref. [4^ . 


Three WS parameterizations [H, are employed 

in our calculations, which were constructed for different 
purposes. The so-called Universal parameter set has been 
extensively used in both nuclear structure (e.g.. Refs. 
[60l - l^ land radioactive decay [63 - 1^ studies. The pa¬ 
rameter set from Ref. (see Table U below and Sec. 
IIIC in Ref. [T^) has been shown to be very successfully 
in both nuclear spectroscopy [68l - [7ll| and binding energy 
[Tsil calculations. In particular, it gives a good descrip¬ 
tion of the nuclear momentum of inertia and high spin 
states. We would also like to emphasize that above pa¬ 
rameter set assumes a rather sophiscated radius of the 
form Ro = ro(l -h 0.116r3(Af - + 0.235 fm. 

The parameter set from Ref. [36l| contains much less 
terms and is primarily fitted to doubly-magic nuclei with 
the restriction that the radius is smaller than 1.3A^/^. 
This potential has not been much tested. We are also in¬ 
terested to see whether it has any major defects in global 
nuclear structure calculations. We label the three param¬ 
eter sets as WSl, WS2, WS3 for simplicity. This should 
not be confused with the WS mass formula as proposed 
in Ref. [ll|. WSl has a very strong isospin dependence 
in SO force with k = kso = 0.86. It is much weaker 
in WS2, where kso is a free parameter with the value 
of 0.16. The parameters of the three WS potentials are 
listed in Table ID 


III. GLOBAL CALCULATIONS 

We have performed systematic calculations for all 
even-even nuclei with Z, N > 8 within the proton and 
neutron driplines as defined in Fig. 12 in Ref. [s^. In 
total 1871 nuclei are calculated and the ground state is 
defined by taking the point with the minimum energy 
within the deformation space {P2,7,134)- We focus in this 
paper on the detailed comparison for known nuclei and 
thus do not explore the positions of the driplines as pre¬ 
dicted by the present model, which may depend on the 
choice of the macroscopic energy model and the uncer¬ 
tainties induced by the pairing interaction The 

calculated potential energy surfaces for all above nuclei 
as well as the deformations and microscopic energies cor¬ 
responding to the global minimum will be available on 
our web page [s^. 


III.l. Systematics on nuclear deformation 

The calculated quadrupole deformations (32 with the 
different WS parameter sets are plotted in Figs. [Tl[5]and 
[31 (32 is assumed to be positive. A prolate (oblate) shape 
corresponds to 7 ^ 0(60) deg. The nucleus would be of 
maximal triaxial deformation if one has 7 ^ 30 deg. The 
7 deformation was not taken into account in some recent 
mass-formula calculations. Its importance is emphasized 
in Ref. (t^. 

In most cases, as can be seen from Fig. |4l the deforma¬ 
tions calculated by the three parameters are quite close 
to each other. For calculations with parameters WSl 
and WS2, there are only 72 cases where differences be¬ 
tween the quadrupole deformation A /32 are larger than 
0.1, among which one has 20 with A /32 > 0.2 and 7 with 
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TABLE I. Parameters of the three Woods-Saxon potentials employed in our calculations. For WSl and WS3 one has Ro — 
whereas for WS2 it is Ro = ro(l + 0 . 116 r 3 (Af — Z)l + 0.235fm and Rso = Ro- 



Vo(MeV) 

ro(fm) 

rso(fm) 

a,aso(fm) 

A 

K 

KSO 

WSl 

49.6 

1.347(n)/1.275(p) 1.31(n)/1.32(p) 

0.7 

35(n)/36(p) 

0.86 

K 

WS2 

53.75 

1.19 

1.19 

0.637 

29.49 

0.791 

0.162 

WS3 

50.92 

1.285 

1.146 

0.691 

24.08 

0.644 

—n 



FIG. 2. (Color online) Same as Fig. [T] but for calculations 
with WS2 from Ref. [m. 
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FIG. 3. (Color online) Same as Fig. [T] but for calculations 
with WS3 from Ref. [M. 


A ^2 > 0.3. The latter correspond to nuclei ®°Sr, 

s^Zr, s^Mo, 292Cn, ^o^IIS, 3O6l20. For the WSl and 
WS3 parameters, there are 118 cases where differences 
between the quadrupole deformation A /32 are larger than 
0.1, among which one has 50 (17) with A /32 > 0.2(0.3). 
As discussed in Ref. [s^ , many of them are neutron-rich 
nuclei with N around 20 and 40 due to the vanishing of 
the shell closure. As can be seen in Fig. [3l for calcula¬ 
tions with WS3, the ground states of a few Sn and Pb 
isotopes are calculated to be deformed, which are slighter 
lower in energy than the spherical minima. 

We also compared our calculated deformations with 
those given by the mass calculations of Refs. ilHl- 
For the deformations given by WS3 parameter and those 
from Ref. Q, there are as many as 66 (34) cases with 
A /32 > 0.2(0.3). However, in the latter cases with 
A /32 > 0.3, 27 are in the superheavy region around 
= 192. In comparison to those given in Ref. [m, 
there are 30 cases where one has A /32 > 0.2. Besides 
the few differences around N = 192, the other cases are 
mainly around N,Z = 14, 28 and 40. In particular, ^®Zr 



FIG. 4. (Golor online) Comparison between the deformations 
calculated from the three different WS parameters. We plot 
the results from parameters WS2,WS3 as a function of those 
from WSl. 



FIG. 5. (Color online) Systematic calculations on the gamma 
deformation 7 of even-even nuclei by using the Woods-Saxon 
potential with WSl. Only results corresponding to nuclei 
with 1/321 > 0.1 are shown. 


is predicted to be largely prolate deformed in our calcu¬ 
lation with WS3 but oblate deformed in Ref. [ll|. ®^Ge 
is predicted to be of triaxial shape with 7 = 41 deg in our 
calculation with WS3. 

The gamma deformations as calculated with the dif¬ 
ferent WS parameter sets are plotted in Figs. [5l[6]and[7l 
Again, the three calculations give practically the same 
results for most nuclei. As can be seen from the fig¬ 
ures, most nuclei are of prolate shape with 7 < 10 deg. 
The nuclei with triaxial and oblate deformations occur 
mostly in regions where transition from spherical to pro¬ 
late shape is expected. In this context, it may be interest¬ 
ing to mention that there is a long history studying the 
origin of dominance of nuclei with prolate shape, which 
may be sensitive to the details of the WS parameteriza- 
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FIG. 6. (Color online) Same as Fig. [5] but for calculations 
with WS2. 



FIG. 7. (Color online) Same as Fig. [5] but for calculations 
with WS3. 

tion [zllzi. However, the predicted deformations seem 
pretty stable within the range of the variety of WS pa- 
rameterizations employed in this work. 

III.2. Deformations of selected isotopes 

As can be inferred from Figs. [TJUl the deformations 
predicted by the three calculations are quite similar in 
most cases. We are particularly interested in the differ¬ 
ences between predictions by the three calculations. In 
Fig. [5] we plotted the predicted deformations of Si, Cr, 
Ge, Zr, Ru and Pb isotopes. 

For Si isotopes, the large difference occurs around 
N = Z = lA. The ground state of ^®Si is calculated 
to be of triaxial shape with (32 = 0.17 and 7 = 22° by 
WSl parameter set. It also gives another two energy min¬ 
ima: the second minimum is calculated to be of nearly 
spherical shape with (32 = 0.07 and the third one is of 
oblate shape with (32 = 0.27. The energy difference be¬ 
tween these three minima is only about 100 keV. It can 
be inferred that the shape of ^®Si predicted by the WSl 
calculation is quite soft. The ground state of this nu¬ 
cleus is calculated to be of oblate shape with (32 = 0.32 
and 0.38 in calculations with WS2 and WS3, respectively. 
A coexisting second minimum with triaxial deformation 
is also predicted in these two calculations. For calcula¬ 
tion with parameters WS2(WS3), the energy difference 
between the two minima is about 300 keV (820 keV). 

For Cr isotopes, the largest difference occurs in 
neutron-rich nuclei around N = 40. The nucleus ®^Cr 


is predicted to be spherical in the calculation with the 
WSl parameter, while it is predicted to be of well de¬ 
formed prolate shape with (32 = 0.27 by the WS3 calcu¬ 
lation. For the calculation with the parameter set WS2, 
the coexistence of spherical and prolate shapes has been 
seen. 

For Ge isotopes, large differences occur in nuclei with 
N>56. But for those nuclei, the ground state shapes 
calculated with the three WS parameter sets are rather 
soft in the 7 direction from —30° to 60°. The (32 values 
are quite close to each other. 

For Zr isotopes, the trends of the evolution of the nu¬ 
clear shape with neutron number calculated with the 
three WS parameter sets are quite similar to each other. 
However, the locations of the transitions are slightly dif¬ 
ferent. It can be seen that the transition from well de¬ 
formed shape to spherical shape around A^ = 40 occurs 
earlier in the calculation with the WS3 parameter than 
in calculations by WSl and WS2 parameters, while the 
evolution from spherical shape to oblate shape around 
N = 56 occurs later in the WS3 calculations as com¬ 
pared to those by WSl and WS2. 

For Ru isotopes, the coexistence of spherical shape and 
oblate deformation appears in nuclei around = 42. 
The isotopes with N = 60 — 66 are predicted to be tri- 
axially deformed with (32 ^ 0.25 and 7 ~ 24°. Then 
transition from triaxial to oblate deformation occurs at 
Af = 68 in calculations by WSl and WS2 parameters 
while triaxial shape remains until A^ = 74 in calculations 
with the WS3 parameter. 

For Pb isotopes, large differences occur in light lead 
nuclei around N = 102 and in heavy nuclei around 
Af = 170. It is well known that the coexistence of spheri¬ 
cal, prolate and oblate shapes occurs in light lead nuclei. 
The ground-state deformations of those nuclei are spher¬ 
ical according to calculations with the WSl and WS2 
parameter sets whereas the prolate minima are lower in 
energy than the spherical minima in calculations with the 
WS3 parameter. This may be due to the simplified treat¬ 
ment of the radius parameters in WS3. For neutron-rich 
nuclei with N ^ 170, the ground-state deformation is cal¬ 
culated to be spherical by WS2 parameter while oblate 
deformation is obtained by the WSl and WS3 calcula¬ 
tions. 

The isospin-dependence of the SO force can indeed lead 
to large fluctuations in ground state deformations for cer¬ 
tain nuclei. Its effect on light neutron-rich nuclei is dis¬ 
cussed in Ref. [1^. To explore this point further, we 
have done schematic calculations with WSl and WS2 by 
changing the sign of the isospin-dependence of the SO 
force Kso- In Figs. l^andlTOl we plotted a few selected 
isotopic chains where the flip is expected to have large 
effects on the shape prediction. It is not surprising to see 
that the effect is larger in calculations with the WSl pa¬ 
rameter which has a strong isospin-dependent SO force 
than that in WS2. However, in the latter case, large de¬ 
viations can also be seen in neutron-rich nuclei including 
those around A^ = 40 and 60. 
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FIG. 8. (Color online) Deformations of Si, Cr,Ge, Zr, Ru and 
Pb isotopes as predicted by different WS calculations. 



FIG. 9. (Golor online) Comparion of ground state deforma¬ 
tions for selected isotopes with and without the flip in nso 
for calculations with WSl. 


III.3. The calculated microscopic energies 


The microscopic energy depends on the deformation 
as well as the single-particle structure of the nucleus to 
be studied. Even though it is not a direct observable, 
the microscopic energy can provide an interesting test 
to our single-particle potentials. In Figs. [TTl [TH and [T 3 ] 
we plotted the total microscopic energy, Edef{N, Z, /3) -|- 
Eshell {N, Z,f3) + Epair {N, Z, f3), for all three calculations 
as a function of N and Z. 



FIG. 10. (Color online) Same as Fig. [9] but for calculations 
with WS2. 
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FIG. 11. (Golor online) Systematic calculations on the to¬ 
tal microscopic energies, Edef{N,Z,l3) -I- EsheiiiN, Z, 13) + 
Epair{N, Z, f3), of even-even nuclei by using the Woods-Saxon 
potential with WSl. 
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FIG. 12. (Golor online) Same as Fig. Illl but for calculations 
with WS2. 
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FIG. 13. (Color online) Same as Fig. 1111 but for calculations 
with WS3. 



FIG. 14. (Color online) Comparison between calculated total 
microscopic energies from the three calculations. We plot the 
results from parameters WS2, WS3 as a function of those 
from WSl. 


For most cases, the three calculations give similar de¬ 
scription on the microscopic energy, as can be seen in 
Fig. [M] For the two calculations with parameters WSl 
and WS3, there are in total 124 cases that the differences 
are larger than 2 MeV. For calculations with parameters 
WS3 and WS2, there are 23 cases with differences larger 
than 2 MeV. For calculations with parameters WSl and 
WS2, there are 28 cases with differences larger than 2 
MeV. 

In some studies, only the shell correction and pairing 
correction energies are considered. The sum of these two 
corrections are plotted in Figs. [131 [H] and [TT] Com¬ 
paring the two calculations with parameters WS3 and 
WSl, there are in total 163 cases that the differences are 
larger than 2 MeV, among which there are 7 cases with 
differences larger than 4 MeV. For calculations with pa¬ 
rameters WS3 and WS2, there are 48 (2) cases with dif¬ 
ferences larger than 2 (4) MeV. Differences are mainly 
seen around magic numbers 20, 28, 40 as well as in su¬ 
perheavy nuclei. For calculations with parameters WSl 
and WS2, there are 66 cases with differences larger than 
2 MeV. 

It should be mentioned that in both cases studied 
above, the energies and deformations correspond to the 
minimum in the total energy surface. 
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FIG. 15. (Color online) Systematic calculations on the shell 
plus pairing correction energy of even-even nuclei by using the 
Woods-Saxon potential with WSl. 
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FIG. 16. (Color online) Same as Fig. 1151 but for calculations 
with WS2. 

III.4. The second minimum 

There is a long history in nuclear physics studying the 
so-called shape coexistence and many prominent exam¬ 
ples have been found [H, [t^, |^. It is also important 
for our study of radioactive alpha [73 and proton decays 
[^ . It is beyond the scope of this paper to analyze in 
detail all cases with possible shape coexistence. In par¬ 
ticular, we are interested to compare the regions where 
shape coexistence is calculated to occur by the three dif¬ 
ferent WS parameter sets. As mentioned above, in a few 
cases the ground state deformations predicted by the dif¬ 
ferent calculations are quite different, which are related 
to the fact that the coexistence of low-lying states with 
different shapes is expected in those nuclei whereas their 



FIG. 17. (Golor online) Same as Fig. 1151 but for calculations 
with WS3. 
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FIG. 18. (Color online) Systematic calculations on the energy 
difference between the second and first deformation minima 
in even-even nuclei by using the Woods-Saxon potential with 
WSl. 



FIG. 19. (Color online) Same as Fig. 1181 but for calculations 
with WS2. 


orders giving by different calculations are different. 

In Figs, [m [19] and |20| we plotted the energy differ¬ 
ence between the second and first minima in the calcu¬ 
lated potential energy surface for even-even nuclei over 
the whole nuclei chart. Only cases with energy difference 
smaller than 1 MeV are selected for simplicity. It can 
be seen from the figures that the three calculations give 
quite similar pattern. Shape coexistence is expected in 
the superheavy nuclei and nuclei around N = 120 and 
170. It may also happen in light neutron-rich nuclei as 
well as nuclei around N, Z = 40. 



FIG. 20. (Color online) Same as Fig. 1181 but for calculations 
with WS3. 


FIG. 21. (Color online) Energy differences between calcula¬ 
tions with and without the Pa and Pa deformations. Those 
with differences smaller than 50 keV are left blank for a clearer 
view. 


III. 5. Effect of octupole and higher-order 
deformations 


Our calculations done above are restricted to the de¬ 
formation space {P 2 , 1 ,Pa) for simplicity. Recent mass 
model calculations were done in the space {P2,I3a,Pq) 
in Ref. [n|, (/ 32 , 7 i/ 34 ) in Ref. [zl and {P 2 : Pz, Pa: P&) 
in Ref. Q. It should be mentioned that there has a 
long quest studying the possible existence of static oc¬ 
tupole and higher order correlations at nuclear ground 
states (see, e.g., a recent experiment [t^ and references 
therein and Ref. [t^ for reviews on earlier works). Both 
the Pa and p(, deformations were taken into account in 
the very successful mac-mic mass model calculations [^. 
Multi-dimensional deformation space calculations have 
also been done for superheavy elements [80l| and in recent 
fission studies . The influence of the Pq deforma¬ 

tion in superheavy nuclei around ^®^No was also studied 
in Ref. [8^. Those correlations may affect the results 
shown above. 

In order to understand the influence of those omitted 
deformation freedoms in our calculations shown above, 
we have done calculations within the axially symmetric 
deformation space {P 2 , Pa, Pa, P 5 , Pe) with the WSl pa¬ 
rameter set. The non-axially symmetric 7 deformation 
is neglected due to computation limitations. In Fig. 
we firstly compared calculations with and without the Pa 
and Pa deformations. It is thus seen that in most cases 
the influence of those two deformations on the total en¬ 
ergy is less than 200 keV. Profound Pa deformation is 
seen in nuclei around N = 134 and Z = 88 and in super¬ 
heavy nuclei around N = 180. This is in agreement with 
earlier calculations Q. 

In Fig. 1211 we compared calculations with and without 
the Pe deformation. The largest influence on total energy 
appears in nuclei around N = 152, in agreement with 
Ref. [S^, and in superheavy elements with N ^ 220. 
The large difference is around 850 keV. 
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FIG. 22. (Color online) Energy differences between calcu¬ 
lations with and without the /Se deformation. Those with 
differences smaller than 50 keV are left blank. 


IV. THE TOTAL BINDING ENERGY 

To explore the description power of the single-particle 
potentials studied above on the (negative) bindiM en¬ 
ergy, we added the liquid drop energy for Emac 
It is chosen to be [l^ 

Eldm = ayA -\- + ^syrn^ (T -h 1) /A 

EasymsT {T+ 1)E + C/^— (13) 

where the terms represent the volume energy, sur¬ 
face energy, symmetry energy, surface symmetry energy. 
Coulomb energy and correction to Coulomb energy due 
to surface diffuseness of charge distribution, respectively. 
The coefficients a„, Og, Ogym, a gyms, C and C 4 are free pa¬ 
rameters to be determined. T = \N — Z\/2 is the isospin. 
As in Refs. [13, the free parameters are firstly de¬ 

termined by minimizing the cr^ value in comparison with 
the experimental binding energies as 

fr" = ^Y.'^Eyaiy{N,Z)-Eyg,pt{N,Z)f , (14) 

* N,Z 

where only nuclei that are heavier than and have an 
experimental error smaller than 100 keV are considered 
and n is the total number of data. The experimental data 
are taken from Ref. [s^. We did not consider the influ¬ 
ence of the uncertainty induced by the errors in experi¬ 
mental binding energies, which is supposed to be small 
as compared with the unknown theoretical uncertainties 

M- 

The parameters thus fitted for the three WS calcu¬ 
lations are given in Table [TTl For calculations with pa¬ 
rameters WSl, one has only 2 cases with deviation be¬ 
tween theory and experimental data larger than 2 MeV. 
Whereas the numbers are 14 and 44 for calculations with 
parameters WS2 and WS3. In all three calculations the 
largest deviations appear at V = 126 isotopes and 
2i6Th. The deviations of the calculated binding ener¬ 
gies from experimental data are plotted in Fig. [53] as a 
function of mass number. All three calculations, in par¬ 
ticular those with parameters WS3, still show some kind 
of systematic deviations around the shell closures, which 


TABLE II. The parameters for the liquid drop model, mean 
deviation and maximum deviation for the different WS cal¬ 
culations as determined by using the least square deviation 
criterion. 



WSl 

WS2 

WS3 

Qjv 

-15.707 ± 0.017 

-15.678 ± 0.017 

-15.704 ± 0.017 

as 

18.302 ± 0.075 

18.197 ± 0.075 

18.450 ± 0.075 

asym 

117.481 ± 0.524 

117.431 ± 0.524 

118.968 ± 0.524 

asyms 

-161.323 ± 2.990 

-161.237 ± 2.990 

-172.024 ± 2.990 

C 

0.717 ± 0.001 

0.716 ± 0.001 

0.719 ± 0.001 

G4 

-0.882 ± 0.072 

-0.896 ± 0.072 

-1.081 ± 0.072 

a 

0.612 

0.731 

1.162 

Max. dev. 

2.320 

2.756 

4.490 


4 


^ 2 

> 

u 

S 

W 0 
<1 

-2 


20 40 60 80 100 120 140 160 180 200 220 240 260 280 

A 

FIG. 23. (Golor online) Deviations of calculations on binding 
energies from experimental data as a function of mass number 
A. A positive value means the binding energy is overestimated 
by theory. 



indicate that the shell correction may have not been fully 
taken into account. 

For binding energies calculated from WS3, in 16 cases 
the deviations from experimental data are larger than 3 
MeV, among which one has 5 cases with N,Z = 20 and 
10 cases with N around 126. For those few nuclei around 
N or Z = 20, the proton and neutron shell gaps pre¬ 
dicted by the WS2 and WS3 parameter sets are much 
larger than those from WSl (by about 2 MeV for neu- 



FIG. 24. (Color online) Same as Fig. [23] but for Pb (dotted 
line) and Po (solid line) isotopes below N = 126 as a function 
of neutron number N. 
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FIG. 25. (Color online) Differences between calculated bind¬ 
ing energies with parameter set WSl and calculations with 
parameters WS2 and WS3 as a function of mass number A. 


tron and 1-1.5 MeV for proton). For the nuclei with 
N ~ 126 and Z between 84 and 92 with large deviations, 
our calculations show that the proton and neutron shell 
gaps predicted by the three calculations are similar to 
each other. However, the neutron shell correction ener¬ 
gies predicted by the WS3 calculation are about 2 MeV 
more attractive than those from WSl. The results from 
WS2 are below those values from WSl and WS3. Ac¬ 
tually, as can be seen from Fig. [Ml the deviations from 
experimental data show a systematic increasing trend as 
a function of neutron number for the different isotopic 
chains with neutron number below JV = 126. The in¬ 
creasing trend is strongest in calculations with WS3 and 
weakest in WSl. For all three calculations, the increas¬ 
ing trend is stronger in the isotopic chains above TV = 82 
than that in Pb isotopes due to the different descrip¬ 
tions of the proton shell corrections. We notice that the 
single-particles energies for levels near the Fermi surface 
given by the WS3 and WSl calculations are very close to 
each other. The difference in the shell correction energy 
may be due to the fact that the single-particle energies 
of the lowest levels are slightly deeper in WS3 than those 
in WSl. The binding energy of ^®Si is underestimated 
by 3.09 MeV in calculations with WS3. This is related 
to the fact that the ground state deformations predicted 
by the WS2 and WS3 calculations are much larger than 
that from WSl, which leads to large repulsive deforma¬ 
tion correction energies -E^ef in those two calculations. 
On the other hand, the shell correction energies given by 
the three calculations are quite similar to each other. 

We are particularly interested in the differences be¬ 
tween different calculations, which may shed light on our 
understanding of the theoretical uncertainty. In Fig. |55| 
we plotted the differences between the three WS calcula¬ 
tions as a function of A for all even-even nuclei considered 
in this work. The largest differences between calculations 
with parameters WS3 and WSl appear in the neutron- 
rich ^°Ca and in superheavy nuclei around Z = 120 and 
N = 254. 

In Figs. [Ml and E71 we compared our calculations on 
the binding energies with mass formula calculations from 



FIG. 26. (Color online) Differences between our mac-mic cal¬ 
culations on binding energies and those of Ref. [Tl| as a func¬ 
tion of mass number A. 



A 


FIG. 27. (Color online) Differences between our mac-mic cal¬ 
culations on binding energies and those of Ref. as a func¬ 
tion of mass number A. 

Refs. sini. In the former case, systematic large devi¬ 
ations are seen in the superheavy region. Whereas the 
deviation of our calculations from those of Ref. 0 seems 
to be much smaller. 

It may also be of interest to determine the parameters 
of the liquid drop model by using the so-called minimax 
fitting procedure. The object of the minimax fit is to find 
the minimum of the maximum deviation (the least worst 
result) as 

£ = argminmax|F;Expt.(^) - Acaic.(A, x)|, (15) 

X A 

where x denote the set of parameters to be determined, 
argmin (argmax) stand for the argument of the mini¬ 
mum (maximum) for which the value of the given ex¬ 
pression attains its minimum (maximum) value within 
a given set of x. The parameters thus determined are 
listed in Table lllll The results given by calculations with 
these parameters are pretty similar to those predicted 
with parameters from Table m This is particularly the 
case for the WSl parameterization. For calculations with 
the WS2 and WS3 parameters, the maximum deviations 
can be largely reduced by applying the minimax fitting 
criterion. However, the corresponding mean deviations 
increase noticeably, particularly for WS3. 
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TABLE III. The parameters for the liquid drop model, mean 
deviation and maximum deviation for the different WS calcu¬ 
lations as determined by using the minimax fitting criterion. 



WSl 

WS2 

WS3 

du 

-15.649 

-15.672 

-15.673 

ds 

17.994 

18.248 

18.266 

dsym 

115.117 

118.397 

116.602 

dsyms 

-148.802 

-172.152 

-160.014 

C 

0.714 

0.7161 

0.715 

C 4 

-0.769 

-0.919 

-0.884 

a 

0.669 

0.957 

1.650 

Max. dev. 

1.980 

2.145 

3.365 


As can be seen from the two tables, most coefficients of 
the liquid drop model can be well constrained by fitting 
to available data except those of the surface symmetry 
energy and surface correction to Coulomb energy. These 
two terms are quite sensitive to the choice of the differ¬ 
ent WS parameterization as well as the different fitting 
criteria. This is expected as they may be more correlated 
with the surface properties of the nuclei induced by the 
nuclear shell structure. 


V. SUMMARY 

We have calculated systematically the microscopic 
energies and nuclear deformations of even-even nuclei 
within the mac-mic framework with three different WS 
parameterizations (denoted as WSl [i^, WS2 and 
WS3 [ 3^ 1 which were constructed primarily for nuclear 
spectroscopy calculations. The first two WS parameteri¬ 
zation have been previously shown to be very successful 
in reproducing many aspects of nuclear structure and de¬ 
cay properties. The simplified WS3 parameter contains 
an unusual SO term and was constructed mainly to ex¬ 
plain the shell evolution in light neutron-rich nuclei . 
We are particularly interested to see whether the isospin 
dependence of the SO force has any global influence on 
the binding energy and on the deformation. 

It is found that the ground state deformations pre¬ 
dicted by the three calculations are quite similar to each 
other. Large differences are seen mainly in neutron-rich 
nuclei and in superheavy nuclei. Systematic calcula¬ 
tions on the shape-coexisting second minima are also pre¬ 


sented. The total binding energy is estimated by adding 
the macroscopic energy given by the usual liquid drop 
model with its parameter fitted by using the least square 
root and minimax fitting criteria. Our calculations on the 
deformation and binding energy are also compared with 
those predicted by two available mac-mic mass formulas. 

It is gratifying to notice that the WS parameters, 
which were fitted only to single-particle states in spher¬ 
ical nuclei, indeed do so well in nuclear mass calcula¬ 
tions. One may speculate that, on a global scale, there 
are smooth effects in the microscopic energy that are 
only depend on the bulk properties of the WS poten¬ 
tial. On the other hand, it is expected that there are 
specific nuclear structure effects that may be sensitive to 
the details of the potential, e.g., the isospin dependence 
of the spin-orbit force and the isospin dependence of the 
radius. While we describe well global properties, it is 
imporant to find the essentials that can differentiate a 
successful microscopic potential from unsuccessful ones. 
In other words, one can identify local properties, e.g., 
the emergence of new subshells and other spectroscopic 
properties, that can help pin down the sign and strength 
of above higher order terms of the WS potential without 
deteriorating its global behavior. 

It is noticed that, in a few cases in heavy nuclei, cal¬ 
culations with the WS3 parameter set can deviate from 
experimental data by as large as 4 MeV (Fig. 23). We 
hope that this can be improved by fine-tuning the pa¬ 
rameterization of the potential (WS3 contains less terms 
than the other two parameter sets employed). We also 
hope that a better understanding of the theoretical un¬ 
certainties of the mac-mic model as well as the WS pa¬ 
rameterization can be obtained. In the future, we are 
also interested to see if one can pin down the isospin de¬ 
pendence of the SO force by fitting to binding energies 
as well as other properties of both stable and unstable 
neutron-rich nuclei. 
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